Random walk on disordered networks 



Tomaso Aste 

(February 1, 2008) 
Equipe de Physique Statistique, LUDFC Universite Louis Pasteur, 
3, rue de I'Univresite, Strasbourg France 
e-mail: tomaso@fresnel.u-strasbg.fr 
PACS: 66.30.-h, 05.40.-j, 61.43. +j 



Abstract 

Random walks are studied on disordered cellular networks in 2-and 3- 
dimensional spaces with arbitrary curvature. The coefficients of the evolution 
equation are calculated in term of the structural properties of the cellular sys- 
tem. The effects of disorder and space-curvature on the diffusion phenomena 
are investigated. In disordered systems the mean square displacement displays 
an enhancement at short time and a lowering at long ones, with respect to 
the ordered case. The asymptotic expression for the diffusion equation on hy- 
perbolic cellular systems relates random walk on curved lattices to hyperbolic 
Brownian motion. 

The simplest disordered networks are space-filling random partitions of space by cells. The 
cells are convex, irregular polygons in two dimension {2d) and irregular polyhedras in 3d. 
Disorder imposes the incidence numbers at their minimum values {d+1 edges incident on a 
vertex in (i-dimensions). These cellular networks are also known in the literature as "froths", 
since the soap-froth is the archetype of such structures. The space tiled by the froth can be 
curved. In this case, the intrinsic dimension of the cellular system (df) does not coincide 
with the dimension (d) of the embedding space. 

Froths are the structures which characterizes a broad class of natural systems such as 
polycrystalline solids, foams, biological tissues etc. Moreover, froths and disordered 

packings are dual systems (see Fig.l). Therefore, amorphous materials, granular solids, 
metallic glasses etc. have structures which are dual of froths 

Many theoretical works, experiments, and computer simulations have been devoted to 
the study of random walk and transport phenomena on disordered P-^j and fractal [^-[15 



systems. Random walk on Euclidean froths is a realistic model for diffusion in disordered 
systems, for signal propagation in granular media and it can be relevant in the study of 
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the evolution properties of natural foams and polycrystalline aggregates. Whereas, random 
walk on hyperbolic or elliptic froths can model transport phenomena in curved spaces. 

Disordered structures are widespread in nature, and natural, disordered, two dimen- 
sional cellular structures can tile curved surfaces (e.g. amphiphilic membranes or epithelial 
tissues). The study of the effect of disorder and space-curvature on the diffusion phenomena 
is therefore of great physical interest and it is the aim of this paper. 

In the present model, the walker starts at time t = from a given cell, then, at each 
(finite) time-step it jumps with equal-probability to one of the neighbouring cells. The radial 
and angular components of the motion respect with the starting cell are decoupled. The 
radial component results the same as of the spherically symmetric model introduced recently 
in |]lT],|l2|. But, in the present case, the diffusion is on realistic cellular systems and all the 
parameters in the evolution equation are given in terms of the properties of the disordered 
cellular structure. 

A d-dimensional disordered cellular froth, can be studied as structured in concentric lay- 
ers of cells at the same topological distance (j) around a given central cell (where the topo- 
logical distance between two cells is the minimum number of {d — l)-dimensional interfaces 
that a path must cross to connect the two cells). The structure is described topologically, 
by two parameters per layer in 2d (number of cells per layer and average coordination in the 
layer), and three parameter per layer in 3d (see [0 for details). 



The number of cells in a layer at distance j from the central cell {K{j)), is related 
to the space-curvature. One finds asymptotically, K{j) oc j'^^'^^, where df is the intrinsic 
dimension of the cellular system. The intrinsic dimension df coincide with the dimension 
d of the embedding space in Euclidean froths (which are cellular tilings of fiat spaces), 
whereas, df > d in the hyperbolic case (tilings of negatively curves spaces) and df < d in 
the elliptic one (tilings of positively curves spaces). A special case, discussed in WE 
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class of hyperbolic froths with K{j) oc exp{(pj). Here the intrinsic dimension diverges. 

Suppose the system being shell-structured-infiatable (SSI) around the central cell. (In 
SSI froths any cell in layer [j) has neighbours in layer (j — 1), (j) and [j + 1) and the layers 
make, concentric, closed rings of cells without "topological defects". See Fig.l and flBl for 
details.) The number of paths connecting different layers can be more easily calculated in 
SSI froths. The extension to the general case of non-SSI froths, follows straightforward. 

Let the central cell be the one where the walker starts at t = 0. Consider at time t the 
walker being in a cell of layer (j) (the cells in layer j are supposed indistinguishable). At 
time t + 1 it has moved outwards to layer (j + 1) or (for j > 0) inwards to layer (j — 1) or 
within the same layer (j), with probabilities PoutU) or Pin{j) or PstayU), respectively. Note 
that, PinU) + Pout{j) + Pstayij) = 1, siucc the walker must move at each time-step. 

The probability P{j,t) of finding the walker in layer (j) at time t is given by 

PU,t)=p,tay{j)PU,t-l) +PoutU - l)P{] - 1) 

+ p,„(j + l)P(j + l,t-l) fori>l, (1) 

and, for j = 0, by 

P(0,t)=p,,„(l)P(l,t-l) . (2) 
With initial conditions P(j, 0) = 5j^Q. 
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The probability Pout{j) is proportional to the number of paths connecting layer (j) with 
layer (j + 1). This number is equal to the number of interfaces (edges in 2d and facets in 3d) 
separating the two layers. Analogously, the probability Pin{j) is proportional to the number 
interfaces between layer (j) and (j — 1). 

In 2d, the number of edges separating layers (j) and (j + 1) is K{j) +K{j + 1) (see fig.l). 
Thus, 



Prn{j)=ji(f)(K{j)+K{j-l) 



(forj>l), (3) 



(and PstayU) = 1 - PoutU) - Pinii) = In Eq.(D we defined, K{0) = and J\f2{j) = 

e{j)K{j), with e(j) the average number of edges per cell in layer (j). For j = 0, one has 
Pout{0) = 1 and pin{0) = 0. 

In 3d the layers are separated by a system of faces which tile a spherical surface: the 
"shell-network" |]T6|]. The number of paths between two successive layers (j) and (j + 1) 



is proportional to the number of facets of the shell-network between these layers. This 



number is 2 '^^-'-'^^|j_^^-' ^ ||16[| , where n(j) is the average number of edges per face in the 



shell-network. We have therefore 



K{j) + K{j + 1) - 8 

' ~ mi) 

PinU) 



PoutU) J\f.^(j) ri(i)-4 



. Kij) + Kij-l)-8 i^orj>l), (4) 



where we defined: K{0) = 2 and Af-s^j) = f{j)K{j), with /(j) the average number of faces 
of the cells in layer (j). For j = 0, we have Pouti^) = 1 and Pm(0) = 0. 

A quantity of interest is the probability 11 that the walker ever returns to the origin. 
This quantity is associated with the mean time spent in the origin (-F(O) = St^o -^(0' ^)) 
the relation 11 = 1 — |TB[. From Eq.(|ip and by using Eqs.(|^) and (^), we obtain 



_ 1 _ 1 

This expression is valid for any froth tiling an unbounded topological manifold. The quan- 
tity J\fd{j)Pout{j) is related to the properties of the structure around the central cell, and 
asymptotically it scales as K{j) (see Eqs.(^) and (D). In a cellular system with intrinsic 
dimension df the number of cells per layer have the asymptotic behaviour K{j) oc 
thus Afd{j)Pout{j) ~ K{j) ^ By substituting into Eq.(^, we obtain II = 1 ioi df < 2, 

and n < 1 for df > 2. This result, already known for random walk on regular lattices, frac- 



tals, trees [0 , and found here in froths, indicates the universality of this critical behaviour 
which is independent of the details of the structure. In Fig. 2 the behaviour of 11 vs. df, 
given by equation (^) for 2d SSI froths with K{j) = K{l)j'^f~^, is shown. 

A quantity generally used to describe the diffusion phenomena is the mean square dis- 
placement (r^)(t) = Y^'jLoj'^P{j,t). The time-dependent diffusion coefficient D{t) is associ- 
ated with this quantity by the relation 2dD(t) = ^{r"^), and the usual diffusion coefficient 
is the limit of D{t) at infinite time. 
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From Eq.(|T]) follows 



+ 1) - (r')(t) = E \Pout{j) + PinU) + 2j [poutU) - P^n{J)] t) . (6) 

When, j 3> 1, and the parameters e(j) = (e), /(j) = (/) and n{j) = (n^) are independent 
of j (this is the expected asymptotic behaviour), Eqs.(0) and (^) give 

for (i = 2 

Pout{j)+Pin{j) = l-Pstay{j) = { = 2Crf , (7) 

^ for = 3 

and, for dj finite, j[Po«t(j) - PinU)] = {df - Thus, from Eq.(§) 

(r2)(t)~2d;Crft . (8) 

The diffusion coefficient is therefore, = ^Cd- Numerical solutions of Eq.(|T]) for 2d and 
3d structures with different intrinsic dimensions and coordinations give diffusion coefficients 
in very good agreement with Eq.(||). Note that, (r^) in Eq.(||) is expressed in term of 
topological distances. The metric quantities can be retrieved from the topological ones by 
multiplying the topological distance j by the average asymptotic distance po between layers. 
For instance, in the hexagonal lattice, po = -^ct, with a the lattice spacing. From Eq.(^), 
one gets therefore (p^) = pKr"^) = a^t, which is the known expression for the mean square 
displacement in the hexagonal lattice. 

The linear dependence on t of (r^) in Eq.(^, indicates a diffusive behaviour. In this case, 
the spectral dimension ds (defined from the exponents (r^) ~ t'^^/^f and P(0, t) ~ t~'^^^'^ ||19|| ) 
coincide with the intrinsic dimension df. 

In disordered froths, topological non-SSI defects are always present. Defects, in layer 



J, are cells which have no neighbours in layer j + 1 (see fig.l and |jT6[). Asymptotically, 



the number of defective cells in layer j is a fraction 5 of the total number of cells K{j). 
Typically, 2d Euclidean, disordered froths have 0.1 < 5 < 0.2 [|l^]. In 2d, the number of 
paths connecting layer j with non-defective cells in layer j + 1 is (1 — 6){K{j) + K{j + 1)) 
(see fig.(l)). Whereas, the number of paths ending in a defective cell is r]6K{j + 1), with 
7] the average number of interfaces added by a defect to the shell between two successive 
layers (typically, 1 < 77 < 1.3 in 2(i Jl^). Therefore, asymptotically, Eqs.(||) and (|^) can be 



extended to the non-SSI case by multiplying expression (0) by the factor (1 - 5 + f ). The 
same result holds in 3d. 

The non-SSI defects have important effects on the froth structure. In particular, it 
has been found in [O that, in 2d non-SSI Euclidean froths, the number of cells per layer 



increases with the distance following a linear law K{j) = Cj + B, with a slope C ~ 9. This 
slope is higher than the both values, C = 27r expected from simple geometrical consideration 
and C = 6 of the SSI hexagonal lattice. An higher increment in the number of cells per 
layer, must correspond to a faster diffusion (higher number of paths to go outward). On 
the other hand, in typical 2d disordered systems, one has {1 — 6 + y) < 1? which indicate 
asymptotically, a slower diffusion in non-SSI froths compared with the ordered SSI case. 
These two opposite behaviours are not in contradiction. Indeed, for j ^ 1 the ratio between 
the number of paths in successive layers depends only on the exponent of K{j) vs. j (i.e. 
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the intrinsic dimension - 1), and not on the slope. Therefore, we expect the diffusion in 
disordered structures compared with ordered SSI lattices, being faster at small distances 
(where the slope of K{j) is relevant) and then becoming slower in the asymptotic part 
(where only the exponent of K{j) is relevant). Fig. 3 shows {r'^)/t calculated numerically 
for a non-SSI 2d Euclidean froth (a) and for the SSI hexagonal lattice (b). The diffusion in 
the disordered structure is faster than in the SSI hexagonal lattice in the first stage (t < 15 
and j < 4), then it slows down to reach an asymptotic behaviour where the mean square 
displacement grows with t more slowly in the disordered than in the ordered system. 

A special behaviour of (r^) is obtained for the 2d SSI hyperbolic froth, studied in Ref. 
| |TB| , which has e(j) = (e) > 6 and K{j) = Cexp{(pj), with (p = cosh~^( ^'^^~^ ). In this case, 
from Eqs.(^, (^) and (^, one derives the asymptotic expression 

(is) -em -2) ^, (9) 

Numerical solutions of Eq.(|I|) for SSI 2d hyperbolic froths with various (e) > 6, give time- 
dependent diffusion coefficient in excellent agreement with expression (|^). The quadratic 
exponent in Eq.(^ indicates a ballistic diffusion and dg = 2df. 

We now write the evolution equation (|l|) in the continuous limit. Let introduce the 
continuous variables p = jpo and r = tro, where po is the average distance between two 
layers and tq is the average time between two jumps. In the asymptotic limit (j = p/po ^ oo 
and t = t/tq ^ oo), when the average topological arrangements of the cells is independent 
of the layer number, equation (|I]) can be written in the continuous form 



{s + 2)K{p)dp 



P{P.r) \ , (10) 



where s is the inflation parameter (s = (e) — 4 in 2d, and s = |((/) — 6)((n^) — 4) — 2 
in 3d [|l6l), which is a quantity associated with the curvature of the manifold tiled by the 
froth (s = 2 Euclidean, s > 2 hyperbolic and s < 2 elliptic froths). Expression ( |ToD is 
the diffusion equation, for a (i- dimensional spherically symmetric cellular system, written 
in polar coordinate. All the information about the structure, its intrinsic dimension, the 
disorder, are contained in the term in the square brackets and in the parameter C^. 

For a random cellular system with finite intrinsic dimension df we have asymptotically 
K{p) ~ p'^-t'~^ and s ^ 2. Therefore the coefficient inside the square brackets in (^) becomes 
and Eq.(|TO|) has the solution 

np,r) = -^^^^^expf--^) ■ (11) 
r(^)(4^C.r)4 V i^QrJ 

The probability P{p,t) increases with p until a maximum at pmax = (2(c?/ — l)^Cdr)^, 
then it decreases exponentially. From solution (0), the mean square displacement results 
(r2) ~ ^P(p, T)dp = 2df^CdT, as already obtained in Eq.(g). 



In a previous paper [|TB[, it was found a class of 2d and 3d hyperbolic SSI froths which 



can be generated iteratively by using a simple recursive equation. These froths have K(p) 



5 



Csinh((y9^), where (y? = cosh" 



is a constant associated with the space curvature (in 



simple 2d cases one can show the equivalence if = v — fc, with k the Gaussian curvature. 
Here, s > 2 and k < 0). For these froths, the term in the square brackets in Eq.(|TDp is 



{s+2) po 



coth(v9^ 



Therefore, the evolution equation ([l0|) takes the form 



d_ 



Pip, 



T 



To dp [ dp 



[s + 2) Po Po 



Pip, 



T 



(12) 



Equation (0) was already known in literature as the diffusion equation in hyperbolic spaces 
with constant, negative curvature 



15,20 



Here the same equation have been obtained 
starting from a tessellation model, making therefore a link between diffusion in curved 
lattices and hyperbolic Brownian motion. At long distances, the coefficient in the square 
brackets in Eg . {\L2^ tends to a constant and the corresponding solution is a Gaussian with 
a constant drift. Here, the probability distribution move ballistically outward with the 
maximum at Pmax = '^f^^dT. This ballistic diffusion is consistent with the result for the 



TO " 



mean square displacement for hyperbolic froth given in Eq.(P). 
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FIGURES 




FIG. 1. A froth is a random partition of space by cells (a). Disorder imposes the incidence 
number at the minimum value (3 edges incident on a vertex in 2d) . Froths are the dual structures 
of disordered packings (b). Such structures can be analysed as organized in concentric layers of 
cells at the same topological distance (j) from a given central cell (j = 0). Some cells (brought 
out by hatcheries in (a)) have neighbours in the internal layer but not in the external one and are 
topological inclusion or "defects"' in the layered-structure. 
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FIG. 2. Probability 11 that the walker ever return to the origin for several values of the intrinsic 
dimension df. The walker always return in the origin when df <2, whereas the probability is less 
than 1 and decreases with df when df > 2. This critical behaviour is independent of the details of 
the structure. (The line is a guide for eyes.) 
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FIG. 3. Mean square displacement (r^) over t vs. time for disordered (a) and ordered (b) 
cellular systems. The average distance of the walker from the starting point is j c± (r^)^/^. At 
short distances (j < 5) the walker diffuses faster in disordered system than in the corresponding 
ordered lattice. Then diffusion in disordered system slows down to reach an asymptotic regime 
where the walker propagates more slowly in disordered system than in the ordered case. 
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